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M. R. Schroeder recently proposed a scheme for compression of motion 
picture data by taking the difference of two successive frames and then 
smearing. 1 The smearing is accomplished by a Hadamard matrix. 

If the Hadamard matrix is of a certain particularly well-understood 
type, then we show that if the input differential picture consists of a small 
odd number of large pulses of identical magnitudes (but arbitrary signs), 
then the output will consist of three components: 

(i) Large pulses of equal magnitude and the correct signs, matching 
each of the input pulses. 

(ii) One additional "stray" large pulse, of magnitude equal to the 
others, but located at a point where the input was zero. 

(Hi) Scattered pulses of amplitude low relative to the pulses of types 
i and ii, but so numerous that they consume (ir — 2)/ir of the total energy 
of the output differential picture. 

We give an explicit formula for the amplitude of each of these pulses. 

The problem of determining the distributions of all possible outputs 
of the proposed system for other classes of inputs is shown to be equivalent 
to the unsolved problem of finding the weight enumerators for the cosets 
of the first order Reed-Muller codes. 

I. INTRODUCTION 

The fact that successive frames of a motion picture are often very 
nearly alike has led to the consideration of schemes which transmit, 
for each point of the picture, the difference between the amplitude 
of the present frame and the amplitude of the previous frame. Since 
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this differential picture* is frequently zero at many points, there is 
reason to hope that the bandwidth required for transmission of the 
differential picture could be greatly reduced by appropriate coding. 
One such coding scheme which has been considered by W. K. Pratt, 
J. Kane, and H. C. Andrews, 2 and refined by M. R. Schroeder 1 is the 
following: let the differential picture be represented by a real n-dimen- 
sional vector, v. (For example, if the picture is represented by a 100 X 
100 grid, then n = 10000.) Let 3C be an n X n Hadamard matrix, 
which is a self-orthogonal real matrix all of whose entries are ±1, 
and let the smeared differential picture (or transformed differential 
picture), x, be defined by x = 3Cv/(n)*. Let Q be the power-preserving 
clipping operator, defined by 

n v - (II x ||)'sgn(x) 
Qx ~ (n)» 

where sgn (x) is the n-dimensional vector whose ith component is 
+ 1 or — 1, depending on the sign of the ith component of x. Since we 
wish the quantizer to have onfy two outputs, we cannot take sgn (0) = 0. 
Unless stated otherwise, we assume that sgn (0) is undefined. Schroeder 
has asserted that the vector y = Qx = Q3Cv/(n)* provides an appropriate 
"encoding" of the differential picture v. To "decode" one computes 
z = l/(n)*3C*y. The question to be studied in this paper is the quality 
of z as an approximation to v. 

Some of the heuristic arguments favoring this proposed scheme are 
the following: Since successive frames are frequently very similar, the 
differential picture will have near-zero amplitude at most points. In a 
typical case when the camera is focused on a moving subject and a 
fixed background, the differential picture will be identically zero at 
all background points. If the subject and the background each has 
uniform color (the simplest plausible case), then the differential picture 
will be nonzero only at those points on the boundary of the subject. 
Furthermore, all of the nonzero amplitudes in the differential picture 
will have equal magnitudes, although their signs will depend on whether 
they are on the leading or trailing edge of the moving subject. The 

* To be precise, the "differential picture" should consist of the difference between 
what the present frame actually is and what the decoder thought the last frame was. 
Since all of the errors in the system are assumed to arise from quantization, rather 
than from any sort of unpredictable noise on the communications channel, the 
encoder may include a replica of the decoder, thereby enabling it to compute what 
the decoder thought the last frame was. Each transmitted differential picture then 
includes an attempt to correct the cumulative effects of all previous errors. In this 
paper, we study only the quantization noise introduced in the encoding and decoding 
of a single differential picture, ignoring the complicated dynamic questions which 
arise when one studies the behavior of the system during several successive frames. 
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conventional manner of encoding the differential picture is to quantize 
the amplitude at each gridpoint. This scheme will introduce no quantiza- 
tion error at all on the background points, which have zero amplitude, 
but a relatively high number of quantization levels may be required 
to keep the quantization errors along the outline of the subject down to 
a tolerable level. The Hadamard transform of the differential picture, 
on the other hand, will have its energy spread out relatively uniformly 
among the grid points. A coarse quantization of the smeared differential 
picture will introduce quantization errors throughout the differential 
picture in a relatively uniform manner. When the quantized smeared 
differential picture is unsmeared, the quantization errors, being some- 
what independent, should tend to cancel out. It is thus hoped that a 
coarser quantization of the smeared differential picture might yield 
a decoded differential picture of the same fidelity as a substantially 
finer quantization of the original, unsmeared differential picture. 

A somewhat more theoretical discussion of the effects of quantization 
in the Hadamard transform domain is given in Section VI of Pratt, 
Kane, and Andrews. 2 The main result is that the Hadamard trans- 
formation preserves energy. Hence, if the amplitudes at the various 
points in the transformed differential picture are independent zero mean 
gaussian random variables, then the energy of the noise introduced by a 
two-level quantizer would be (tt - 2)/x of the total energy in the output 
differential picture, both before and after unsmearing. From this view- 
point, the major attraction of smearing is that it distributes the quantiza- 
tion noise uniformly throughout the picture. If our simple model of a 
differential picture (which lias nonzero amplitude only along the 
outline of the subject) is coarsely quantized, then all of the quantiza- 
tion noise appears on the outline, where it will tend to blur the subject. 
However, if this differential picture is smeared, coarsely quantized, 
and unsmeared, then its quantization noise should be evenly distributed 
throughout the subject and the background. 

We shall now study the relationship between the original differential 
picture, v, and the decoded differential picture, z. It is clear that the 
energy in the vector z is always identical to the energy in the vector v. 
Hence, for purposes of analysis, it is easiest to compute z according 
to the formula 

z = AX.' sgn 3Cv 

where for each frame A is a non-negative scalar chosen to make the 
energy in z equal to the energy in v. In this paper, we often omit the 
actual calculation of A. 
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In the case where v has only one nonzero component, then sgn 3Cv = 
OCv and z = Aac'aCv. Since 3C'3C = nl (where I is the n X n identity 
matrix) it follows that z = v. In other words, the system transmits 
a single pulse without error. 

On the other hand, when v has only two nonzero components, then 
the component with the larger amplitude dominates the component 
with the smaller amplitude. In this case z again has a single nonzero 
component, even though v had two nonzero components. However, 
the choice v = [1, 1 — e, 0, 0, 0, • • • , 0] results in an ambiguity. If 
we instead write v = [1, 1, e, 0, 0, 0, • • • , 0], then we may actually 
find that, in the limit as e -> 0, z -> [1, 1, 0, 0, ■ • • , 0]. If v = [1, 1, 0, 0, 
• • • ,0], then z is undefined because it depends on sgn 0, which is either 
plus or minus. In fact, z is undefined whenever v has an even number 
of nonzero components, all of equal magnitude; but this difficulty might 
be removable by adding an appropriate background noise function 
into v, or by choosing sgn(0)'s independently at random. To avoid 
the necessity of such considerations, we devote our primary attention 
in this paper to the case in which v has an odd number of nonzero com- 
ponents, all of unit magnitude (but arbitrary sign). In this case, every 
component of 3Cv is an odd integer. Since every component of 3Cv 
must therefore have magnitude at least 1, the sgn function is defined 
and the analysis remains valid in the presence of a small background 
noise in any or all components of v. 

II. HADAMARD MATRICES 

The requirement that any three rows of a Hadamard matrix be pair- 
wise orthogonal leads to the immediate conclusion that if n > 2, then 
an n X n Hadamard matrix can exist only if n is a multiple of 4. The 
question of whether or not there actually do exist n X n Hadamard 
matrices for all n = mod 4 is now one of the most intriguing unsolved 
problems in combinatorial theory. Many ingenious constructions have 
been proposed, and several of them succeed in obtaining Hadamard 
matrices for an infinite number of (scattered) values of n. For example, 
if n is a multiple of 4 and n — 1 is a prime-power, then a well-known 
construction based on quadratic residues in the finite field GF(n — 1) 
yields an n X n Hadamard matrix. Many other constructions for 
Hadamard matrices are given in Chapter 14 of Hall, 3 and more recent 
constructions have been presented by Spence, 4 Goethals and Seidel, 
and Wallis. 0,7 The smallest value of n = mod 4 for which no n X n 
Hadamard matrix has yet been constructed is n = 188. 
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For many values of n, there exist Hadamard matrices with additional 
structure. For example, some Hadamard matrices have the property 
that the first row and first column consist entirely of +l's, and the 
remaining (w — 1) X (n — 1) submatrix has the property that each 
of its rows is a cyclic shift of the previous row. Such matrices are called 
cyclic Hadamard matrices. They are known to exist whenever n — 1 
is prime, or when n — 1 is the product of twin primes, or when n is a 
power of 2. A computer search by Thoene & Golomb 8 and some cal- 
culations by Baumert 9 have shown that no cyclic Hadamard matrices 
of other orders less than 1000 exist, with the possible exceptions of 
n = 400, 496, 628, 652, 784, 976. 

From the viewpoint of an algebraic coding theorist, a shortened 
Hadamard matrix (obtained from a standard Hadamard matrix by 
multiplying each row by an appropriate sign to make the first column 
all -f-l's, and then deleting the first column) is equivalent to an equi- 
distant binary code. The n codewords are taken as the rows of the 
shortened Hadamard matrix, with each +1 replaced by and each 
— 1 replaced by 1. Since the dot product of any pair of rows in the 
shortened Hadamard matrix is -1, the distance between any pair of 
words in the binary code is (n + l)/2. Further discussion of such codes 
is given in Section 13.5 of Berlekamp. 10 

The best-understood class of equidistant binary codes is the maximal- 
length shift-register codes, which are also called shortened first-order 
Reed-Muller codes. In addition to being cyclic and equidistant, they 
are linear over the binary field, which means that the component by 
component binary sum of any pair of codewords is another codeword. 
Stated in terms of the original Hadamard matrices, this property means 
that the componentwise product of any pair of rows of the Hadamard 
matrix is another row of the same Hadamard matrix. Although Hada- 
mard matrices with this property are relatively rare, they exist for 
every n which is a power of 2. Because of their correspondence to the 
Reed-Muller codes, these matrices are comparatively well-understood, 
and we shall henceforth confine our discussion to Hadamard matrices 
of this type. Such matrices may be taken as cyclic. 

III. THE INDUCED COORDINATIZATION 

A 2* X 2* Hadamard matrix corresponding to a Reed-Muller code 
induces a (non-unique) coordinatization on the 2 A components of each 
row, associating each component with a /e-dimensional vector over 
GF(2). A set of 2' coordinates is said to form an affine subspaceiff the 
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corresponding 2' fc-dimensional vectors form an affine ;'-dimensional 
subspace over GF(2). Similarly, a set of m coordinates are said to be 
linearly (or affine) independent iff the corresponding fc-dimensional 
vectors are linearly (or affine) independent. 

If a set of m binary fc-tuples a, , a 2 , • • • a m are affine independent, 
then they span an (m — 1) -dimensional affine subspace, each element 
of which has a unique representation of the form 

2/ + 1 

for some j. A set of binary vectors are affine independent iff no subset 
containing an even number of vectors sums to zero. An affine basis for 
the set of all 2* binary fc-dimensional vectors may be selected in various 
ways. The "standard" basis consists of the all-zero vector and each 
of the A: "unit" vectors. In general, any k + 1 affine independent vectors 
may be chosen as a basis. 

If fc is very large, then the probability that a randomly chosen set 
of k + 1 fc-dimensional binary vectors will be affine independent is 
fj« =i (1 _ 2"*), or about 29 percent. The probability that k randomly 
chosen vectors will be affine independent is about 58 percent; for 
k — 1, it is 76 percent. If m « k, then almost every set of m different 
fc-dimensional binary vectors is affine independent. 

The first row of the 2* X 2 fc Hadamard matrix may be taken as all 
+ l's. The 2*" 1 +l's in each of the other rows occur in the components 
corresponding to some (fc — l)-dimensional subspace of the fc-dimen- 
sional binary vectors, and the — l's occur in the components correspond- 
ing to the complementary (fc — l)-dimensional affine subspace. 

The coordinatization induced by the Hadamard matrix is invariant 
under all changes of affine basis. When translated into coding ter- 
minology, this is equivalent to Theorem 15.35 of Berlekamp. 10 

IV. MAIN RESULT AND DISCUSSION 

4.1 Theorem 

Let v be a 2 fc -dimensional vector whose only nonzero components 
are 2m + 1 units occurring at components corresponding to fc-dimen- 
sional binary vectors a n , a x , a 2 , ■ ■ ■ a 2m which are affine independent. 
Let the vector z be denned by the equation 

z = 3C' sgn 3Cv. 

Then the value of z p , the component of z corresponding to the fc- 
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dimensional binary vector /3, is given by 

[ if /3 is not in the affine subspace spanned by 



2fl = i 



«0 , Ct\ , Oi 2 



S2 k - 2m (2m - 2j)\ (2j)\ . 



j] ml (m — j)\ 
and the sign is given by 



if /»= !>,, 



2/ + 1 



4.2 Remarks 

We first notice that the answers do not significantly depend on A;, 
but on 2m + 1, the number of units in the input vector v. 

Since z is defined by an equation of the form z = 3C'y> where the 
energy of y is n, the energy of z is n~ — 4 . For i = 0, 1, ■ ■ ■ 2m we have 



2,. • " = 



H 2 ;)]'- 



For large w, | z ai | 2 is closely approximated by 2 2k /wm. Thus only 
2/t of the total energy in z is located in those components in which v 
has units; the remaining {it — 2)/ir of the total energy is distributed 
throughout the affine subspace. For example, if m = 5, h — 10, wc 
have the following output: 



Input Value | 




Corresponding 
Output Value 


/OutputV 
VValue / 


Number of such 
Coordinates 


1 







252 
28 
12 

252 


16 X 3969 
16 X 49 
16 X 9 
16 X 3969 


(V) = ii 

(V) + (V) = 220 
(V) + (V) = 792 

SI) = i 


Totals 




16 X 2 16 


1024 



For all m, we notice that if /3 = 2^<-o a t > then 



In other words, if the input to the coding-transmission-decoding system 
consists of 2m + 1 pulses of equal magnitudes (and arbitrary signs) 
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located at positions which are affine independent then the output 
may be written as the sum of the following three terms: 

(i) 2m + 1 pulses of equal magnitude (and correct signs) matching 
the input. 

(ii) One stray pulse of the same magnitude as the 2m + 1 correct 
ones. 

(in) Other errors scattered throughout the affine subspace, having 
maximum amplitude [l/(2m - 1)] times as large as the correct pulses. 

Since each of the errors of type Hi has a relative amplitude approaching 
zero for sufficiently large m, one might consider the proposed system 
"successful" in some sufficiently broad sense of that term even though 
these errors consume (v — 2)/ir of the energy in the output signal. 
The error of type ii poses a different problem, even though it consumes 
a negligible fraction of the energy. Further research may be required 
to decide whether these difficulties might be removed by replacing the 
operator Q by another quantizer with more levels. Further study will 
also be required to determine how these quantization errors propagate 
in successive frames in a dynamic system. (See footnote on page 970.) 

V. RELATIONSHIP TO THE WEIGHT ENUMERATION PROBLEM FOR RM COSETS 

In the previous sections we calculated the vector z = 3C'Q3Cv for 
certain specific choices of v. These vectors v were chosen to be "typical" 
in some intuitive sense, and yet sufficiently simple in form to enable 
us to carry through the calculation in closed form, even when the 
dimensions of the 3C matrix (n X n) were large. 

Instead of assuming some ad hoc form for the vector v, we might 
instead ask, what is the range of the operator 3C'Q3C? Except for the 
scalar factor, this is equivalent to determining the range of the operator 
3C' sgn. For, if there exists a vector x such that z = 3C ( sgn x, then 
n sgn x = sgn 3Cz and z = (3C C sgn Xz)/n. Hence, every vector in the 
range of 3C,' sgn is proportional to a vector in the range of JC'QOC, and 
every vector in the range of 0C'Q3C is a stationary point of this operator. 
Stated another way, (3C'Q0C) 2 = 3C'Q3C. 

In more practical terms, an investigation of the vectors in the range 
of 3d' sgn is actually an investigation of the ensemble of possible dif- 
ferential pictures which the proposed system might produce as output. 
This set is identical to the set of differential pictures which the system 
will encode and decode with zero error. 
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If 3C is an ft X n Hadamard matrix, then there are 2" vectors in the 
range of 3C C sgn. For reasonable values of n, 2 n is so large that a complete 
listing of all of these vectors is not feasible. Fortunately, however, 
these 2" vectors fall into a relatively small number of classes, each class 
consisting of those vectors which have the same distribution of mag- 
nitudes of component amplitudes. The problem of determining the 
possible distributions of magnitudes of the component amplitudes of 
a vector in the range of 3C* sgn turns out to be identical to the problem 
of determining the weight enumerators for the cosets of the Reed- 
Mull er code. This equivalence is seen as follows: If y is a real vector 
whose components have unit magnitude, then the number of components 
of 3C'y with magnitude | A | is the number of rows of 3C whose dot 
product with y is ±A. On the other hand, if we convert 1 to and 
-1 to 1, changing 3C to 9 and y to R, then the weight of the binary 
vector sum of R and each codeword of the extended max-length feed- 
back shift register code gives the distance between the received word 
R and the corresponding codeword, and the enumeration of all of these 
weights for a particular R is the weight enumerator for the coset con- 
taining R. If c is a binary codeword for which w(c + R) = d and 
if v is the real vector of dbl's corresponding to c, then y and y disagree 
in d components and agree in n — d components. Therefore, 

y y = n — 2d. 

Since the first order Reed-Muller code contains both the codewords 
in the extended maximum length feedback shift register code and their 
complements, there is a one to one correspondence between RM cosets 
with weight enumerator d , d, , • • • , d n and real vectors in the range 
of 3Q. 1 sgn having magnitudes of component amplitudes distributed as 
follows: di components with amplitudes ±(n — 2i) for i = 0, 1, 2, 
• • • , n/2 — 1, and d n/2 /2 components with amplitude zero. 

The coset weight enumerators for first order RM codes of lengths 
up to 16 have been determined by R. Dick and N. J. A. Sloane. 11 
Their results, and the corresponding distributions of magnitudes of 
amplitude components of the output of the differential picture encoding- 
decoding system are shown in Table I. Those rows which are predicted 
by our main theorem have been checked, and the relevant values of 
in have been listed. 

The coset weight enumerator for the first order RM code of length 
32 was determined by Bcrlekamp and Welch, 12 and the results are shown 
in Table II. 

The coset weight enumerators for first order RM codes of length 
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^64 have not yet been determined. This problem definitely merits 
further research. 
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APPENDIX A 

A Sketched Proof of Main Result 

One of the implications of our main theorem is that if /3 lies outside 
of the affine subspace spanned by the relevant a's, then z fi = 0. A 
generalization of this result is the following: 

Theorem: If the only nonzero components of v all lie in a (k — 1) 
dimensional affine subspace and (3 lies outside of this subspace, then z = 0. 

Proof: Using the elementary properties of RM codes and affine 
subspaces, the original 2* X 2* Hadamard matrix may be partitioned as 



3C = 






9 
L9 -9. 

where 9 is a 2*" 1 X 2 fc_1 Hadamard matrix. In terms of this partition, 
the last 2* -1 components of v are zero and 3Cv is of the form [w, w]' 
for some appropriate 2*~ '-dimensional vector w. We then obtain 

X 1 sgn [w, w]' = 2[z, 0]' 

where 

z' = 9' sgn w'. 

By repeated application of this theorem, we deduce that z$ = 
unless /3 is in the affine subspace spanned by the coordinates of the 
nonzero components of v. If the coordinates of the nonzero components 
of v span a d-dimensional affine subspace, then an appropriate change 
of coordinates allows us to work with a 2 d X 2 rf Hadamard submatrix, 
which is also Hadamard. The original output z vector merely gains a 
factor of two for each omitted dimension. 

Applying these arguments to the main theorem of the text allows us 
to confine our attention to the case when k = 2m. 

Since the RM code is invariant under the full affine group, there is 
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Table II — Coset Weight Enumerators for First Order 
RM Codes of Length 32 











Value 


Boolean 








of m 


function 


Niunber of 






if pre- 


for coset* 


such cosets 


Weights 




dicted 






2 4 6 8 10 12 14 


16 




Even Cosets 




32 30 28 26 24 22 20 18 (halved) 




2345 


496 X 1 


1 15 


16 




2345&12 


496 X 120 


0000220 14 


14 




2345&23 


496 X 35 


0001030 12 


16 




2345&23&45 


496 X 28 


0000060 10 


16 


2 


2345&12&34 


496 X 840 


0000028 14 


8 




2345&123 


17360 X 2 


0010003 16 


12 




2345&123&12 


17360 X 24 


1 1 4 14 


12 




2345&123&24 


17360 X 18 


0000204 16 


10 




2345&123&14 


17360 X 192 


0000124 14 


11 




2345&123&45 


17360 X 32 


0000044 12 


12 




2345&123&12&34 


17360 X 72 


0000044 12 


12 




2345&123&14&35 


17360 X 576 


0000028 14 


8 




2345&123&12&45 


17360 X 96 


1 8 16 


7 




2345&123&24&35 


17360 X 12 


12 16 


4 




2345&123&145 


13888 X 320 


1 1 6 15 


9 




2345&123&145&45 


13888 X 32 


0001006 15 


10 




2345&123&145&24&45 


13888 X 480 


0000036 13 


10 




2345&123&145&35&24 


13888 X 192 


1 10 15 


6 




123 


155 X 8 


10 7 


24 




123&45 


155 X 512 


0000040 28 







123&14 


155 X 168 


00002080 


22 




123&14&25 


155 X 336 


000000 16 


16 




123&145 


868 X 32 


1 1 30 







123&145&23 


868 X 320 


000010 12 


19 




123&145&24 


868 X 480 


0000040 28 







123&145&23&24&35 


868 X 192 


000000 16 


16 




12 


1 X 155 


00004000 


28 


1 


12&34 


1 X 868 


000000 16 


16 




— 


1 X 1 


10 


31 






* These functions are written in an abbreviated notation. For example, the 
second line, 2345&12 indicates that this equivalence class of cosets includes one 
coset whose members are the 64 Boolean functions of the form X2X3X4X5 + 
X1X2 + AX. + BX 2 + CX 3 + DX 4 + EX B + F, where A, B, C, D, E, and F 
are arbitrary binary elements. 



no loss of generality in assuming that a = 0, that a t , a 2 , • • • , a 2m 
are unit vectors, and that 

V ai = +1 for I = 0, 1, ••• , 2m. 

Any other case can be reduced to this case by an appropriate affine 
transformation of coordinates. 

We now determine the distribution of the vector 

x = 3Cv. 
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Table II— Cont'd 







1 3 5 7 9 11 13 


15 




Odd Cosets* 




31 29 27 25 23 21 19 


17 






32 X 1 


10 


31 




12 


32 X 155 


13 


28 




12&34 


32 X 868 


6 10 


16 




123 


4960 X 1 


10 7 


24 




123&12 


4960 X 7 


10 3 4 


24 




123&14 


4960 X 84 


112 6 


22 




123&45 


4960 X 64 


3 17 


21 




123&14&25 


4960 X 336 


6 10 


16 




123&12&45 


4960 X 448 


1 3 13 


15 




123&12&34 


4960 X 84 


2 4 4 


22 




123&145 


27776 X 10 


1 1 12 


18 




123&145&12 


27776 X 6 


1 1 10 


20 




123&145&23 


27776 X 80 


10 3 9 


19 




123&145&45&23 


27776 X 16 


1 1 15 


15 




123&145&24 


27776 X 180 


2 2 10 


18 




123&145&24&23 


27776 X 240 


15 7 


19 




123&145&35&24 


27776 X 240 


1 3 13 


15 




123&145&35&24&23 


27776 X 192 


6 10 


16 




123&145&45&35&24&23 


27776 X 60 


4 16 


12 





f All functions representing odd cosets also contain the term 12345, which is not 
shown in this table. 



Since x is the sum of the all-ones vector (corresponding to the column 
of 3C associated with a = 0) and 2m columns of 3C which correspond to 
linearly independent codewords in the RM code, it is readily seen that 
there are ( 2 ( m ) components of x which have the value 2m + 1 — 2/. 
It follows that y = sgn x has XXo CD +l'a and £E«*i (?) -l's. 

For convenience, we may partition the components of y into 2m -f- 1 
subsets, each of which corresponds to the components of x with the 
same value. We call the set which consists of ( 2 ™) components where x 
had value +1 the "significant" set of components. The sets of in- 
significant components may be matched up in pairs; the set where x 
had value 2m + 1 — 2t being matched with the set where x had value 
2m -f- 1 + 2t. Each of these matched sets contains ( 2 , m ) components. 

We now let /3 be a typical 2m-dimensional binary vector, which is the 
sum of 2j+l as. By an appropriate permutation of basis vectors, 
we may assume that 

2j 2j 

i - 1 i=l 

The equation z = 3C*y now allows us to compute z & as the dot product 
of a particular row of X* and y. The equation /3 = S-L, a { and the 
correspondence to RM codes allows us to express the particular row 
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of 3C' under consideration as the componentwise product of the first 
2j rows of 3C*, ignoring the zeroth row, which is all plus one. Symbolically, 
if we let r x , r 2 , • • • , r 2m denote the first 2m rows of 3C', then 

H = ( r i <8> r 2 ® ' • ■ ® r 2j )-sgn (l>,J 

where "(x)" denotes the componentwise product. Since 2j is even, we 
may also write 

H = K-rO ® (-r 2 ) <g> • • • (g) (-r 2; -)]-sgn (£r,)- 

The dot product is the sum over all 2 2 '" components of the component- 
wise product of (rj (x) r 2 (x) • • • (x) r 2/ ) and y = sgn (S?-i r <)- Since 
there is cancellation of the summands coming from matched sets of 
components into which we partitioned y, we need only consider the 
Cm) "significant" components. On each of these, y takes value +1, and 
the problem reduces to the following: Given a 2m X Cm) matrix, whose 
columns represent all ways of distributing m plus ones and m minus 
ones among 2m rows, compute the sum of the entries in the component- 
wise product of the first 2j rows of this matrix. The solution is obtained 
by noting that if there are i minus ones in the first 2j rows, then the 
componentwise product is (—1)* and this happens in Ci)Cm-V) columns. 
Therefore, 



= ?(-^( 2 ::?)- 



Having already explained the other factors in the more general version 
of the theorem stated in the text, the theorem is reduced to the identity, 

V r-lv(W 2m ~ 2j ) = (-l)'(2m-2j)!(2j)! 

Y { ' \il\ m- i I j\ m! (m - j)\ 

Multiplying through by (ml) 2 /(2m — 2j)l(2j)\ reduces this to the 
equivalent identity, 



^M% m -^ l ^iD 



ij - %i \j 

whose proof is given bj r Riordan (p. 14, line 7 from bottom). 13 Q.E.D. 



APPENDIX B 



An Example 

Suppose that the differential picture consists of a 4 X 4 grid, the 
points of which are lettered as follows: 



REDUCING -MOTION PICTURE BANDWIDTH 
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A 


B 


C 


D 


E 


F 


G 


H 


I 


J 


K 


L 


M 


N 





P 



The signs of the units in the 16 X 16 cyclic Hadamard matrix with 
which the differential picture is smeared may be taken as: 



0C = 



4- 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


— 


+ 


+ 


- 


— 


+ 


— 


+ 


- 


— 
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+ 
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+ 


— 


- 


+ 


— 


+ 


— 


- 


— 


— 


+ 


+ 


+ 


— 


+ 


+ 


— 


- 


+ 


- 


+ 


— 


— 


— 


— 


+ 


+ 


+ 


- 


+ 


+ 


— 


— 


+ 


— 


+ 


- 


- 


- 


- 


+ 


+ 


+ 


+ 


+ 


+ 


- 


— 


+ 


— 


+ 


- 


- 


— 


- 


+ 


+ 


+ 


- 


+ 


+ 


— 


— 


+ 


— 


+ 


- 


— 


- 


- 


+ 


+ 


+ 


- 


+ 


+ 


- 


— 


+ 


— 


+ 


- 


- 


- 


— 


+ 


+ 


+ 


- 


+ 


+ 


+ 


- 


+ 


- 


+ 


- 


- 


- 


- 


+ 


+ 


+ 


- 


+ 


+ 


- 


+ 


+ 


— 


+ 


— 


- 


- 


— 


+ 


+ 


+ 


— 


+ 


+ 


- 


— 


+ 


— 


+ 


— 


- 


- 


- 


+ 


+ 


+ 


- 


+ 


+ 


- 


- 


+ 


+ 


+ 


— 


- 


- 


- 


+ 


+ 


+ 


— 


+ 


+ 


- 


- 


+ 


— 


+ 


- 


— 


- 


- 


+ 


+ 


+ 


— 


+ 


+ 


- 


— 


+ 


- 


+ 


+ 


- 


— 


— 


+ 


+ 


+ 


- 


+ 


+ 


- 


- 


+ 


- 


+ 


- 


+ 


- 


— 


+ 


+ 


+ 


- 


+ 


+ 
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- 


+ 


- 


- 


- 



The induced coordinatization may be read from the 2nd, 3rd, 4th and 
5th rows. It is: 
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Representation as Sum of 






Odd Number of Affine 


Grid Point 


Coordinatization 


Basis Vectors 


A 





H + I + L 


B 


1 


E+H+I+L+P 


C 


10 


H + L + P 


D 


10 


E + I + L 


E 


10 1 


E 


F 


11 


E + H + L 


G 


110 


E + L + P 


H 


110 1 


H 


I 


10 10 


I 


J 


10 1 


I + L + P 


K 


10 11 


E + I + P 


L 


111 


L 


M 


1111 


H + I + P 


N 


1110 


E + H + I 





110 


E + H + P 


P 


10 


P 



The coordinates of B through P may also be taken as the successive 
contents of this shift register: 



Now suppose the differential picture is this: 















+1 








-1 


+1 








-1 











-1 
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The coordinates of the nonzero inputs are as follows: 

E 10 1 

H 1 1 1 

I 10 10 

L 1 1 1 

P 10 

1= sum 

These are affine independent, so our main theorem applies. The stray 
output pulse will be located at point B, since this is the point whose 
coordinates are the vector sum of the coordinates of the other inputs. 
We now calculate the output vector step by step, without using the 
theorem. The input picture is: 

v = [0, 0, 0, 0, +1, 0, 0, -1, +1, 0, 0, -1, 0, 0, 0, -1]. 

The "smeared picture" is x = 3Cv, which is given by the sum of the 
following rows: 

+ - + + -- + - + ---- + + + 

- + + - + - + + + + --- + -- 
+ - + - + ---- + + + - + + " 
-- + + + + --- + -- + + - + 

- + --- + -- + + - + - + + + 



x = -l,-l,+3, -1, + 1,-1,-1, -3, + l,+3,-3,-l,-3,+5, + l, + l 
The quantized smeared picture is 
y = Qx 

= [ + - + + + + + +]• 

The decoded differential picture is given by 

z = 3Cy 

= [-2, -6, +2, +2, +6, -2, -2, -G, +G, -2, +2, -G, -2, 
+ 2, -2, -6]. 
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When scaled down by a factor of six and placed on the grid, the output is 



_1 

3 


-1 


+ ^ 


+ * 


+ 1 


1 
3 


1 
3 


-1 


+ 1 


1 
3 


+* 


-1 


1 
3 


+ 4 


1 
3 


-1 



In order to match the total power of the input, the output must be 
scaled down slightly more. 
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